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Abstract. We present results of the analysis of data collected in 57-ks XMM-Newton and 505-ks Chandra obser¬ 
vations of the near by (~ 120 pc) isolated neutron star RX J1856.5-3754. We confirm most of the statements made 
by Burwitz et al. (2001) who discussed the original 55-ks Chandra data. Detailed spectral analysis of the com¬ 
bined X-ray and optical data rules out the currently available nonmagnetic light and heavy element neutron star 
atmosphere (LTE) models with hydrogen, helium, iron and solar compositions. We find that strongly magnetized 
atmosphere models also are unable to represent the data. The X-ray and optical data show no spectral features 
and are best fitted with a two-component blackbody model with kT{^ x — 63.5 eV and x — 4.4 (d/120 pc) km 
for the hot X-ray emitting region, and fcTi/g < 33 eV and > 17 (d/120 pc) km for the rest of the neutron 

star surface responsible for the optical flux. The large number of counts collected with XMM-Newton allows us 
to reduce the upper limit on periodic variation in the X-ray range down to 1.3% (at a 2a confidence level) in 
the 10“® — 50 Hz frequency range. In an attempt to explain this small variability, we discuss an one-component 
model with kT^ ~ 63 eV and ~ 12.3 (d/120 pc) km. This model requires a low radiative efficiency in the X-ray 
domain, which may be expected if the neutron star has a condensed matter surface. 
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1. Introduction 

1.1. Thermal radiation from isolated neutron stars 

Thermal emission from the surface of an isolated neutron 
star (NS) could be very useful for determining the star’s 
mass M and the radius R, what in turn puts important 
contraints on the equation of state of matter at supranu¬ 
clear densities. 

Specifically, the detection of spectral features in radi¬ 
ation of isolated NSs may provide: 


RXJ1856.5-3754 (or RXJ1856 throughout the text) is the 
brightest of them in X-rays. 

We use apparent temperatures and radii 
which are measured by a distant observer, throughout 
this paper. The true parameters as measured at the neu¬ 
tron star surface are given by '?bb = T^ (1 — rg/i?)“^/^ 
and i?bb = Rbb(l “ rglRyi"^ where rg = 2GM/c^ is the 
Schwarzschild radius of the NS. 

1.2. Previous Observations 


— chemical composition of the NS surface, 

— the surface gravity from the line broadening, 

— the mass-to-radius ratio from measurement of gravita¬ 
tional redshift of the lines, 

— strength of the surface magnetic field. 


ROSAT was the first satellite with a sufficient sensi¬ 
tivity in the soft X-ray band (0.1 — 2.0 keV) to start a 
systematic search and study of isolated NSs which reveal 
blackbody-like thermal emission. Seven such objects emit¬ 
ting soft thermal X-ray spectra, having high X-ray/optical 
flux and invisible in radio band were detected with ROSAT 
(see, e. g., Zampieri et al. 2001 for the complete list). 


Send offprint requests to: Vadim Burwitz, 
e-mail: burwitz@mpe.mpg.de 


Walter et al. (1996) and Neuhauser et al. (1997) first con¬ 
cluded that RXJ1856 is an isolated NS. The ROSAT data 
showed that RXJ1856 has a non-variable X-ray flux of 
~ 1.5 X 10“^^ erg s“^ cm“^ and a soft spectrum with a 
blackbody temperature of kT^ ~ 57eV. Very interest¬ 
ingly, the Planckian shape of the spectrum was confirmed 
later with the fi rst hi gh-resolution Chandra LETGS data 
(Burwitz et al. 2001 ; hereafter BOl). Actually, the grat¬ 
ing spectrum did not show any significant deviations from 
a Planckian, which are expected in spectra of NS atmo¬ 
spheres containing heavy chemical elements. 

The faint optical counterpart with F ~ 26 mag of 
RXJ1856 was discovered by Walter & Matthews (1997) 
and confirmed by Neuhauser et al. (1998). Walter (2001) 
and Neuhauser (2001) used the faintness of the opti- 
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Table 1. Journal of Chandra and XMM-Newton observations 


Obs. Id. 

# 

Instrument (mode) 

Obs. Start 
[UT] 

Obs. End Exposure 
[UT] [ksec] 

Obs. 

Type 

Remark 

Chandra 

113 

LETGS 

2000/03/10 

2000/03/10 

55.48 

GTO 


3382 

LETGS 

2001/10/08 

2001/10/09 

101.96 

DDT 


3380 

LETGS 

2001/10/10 

2001/10/12 

167.50 

DDT 


3381 

LETGS 

2001/10/12 

2001/10/14 

171.13 

DDT 


3399 

LETGS 

2001/10/15 

2001/10/15 

9.32 

DDT 




total exposure 

505.39 




XMM-Newton 






0106260101 

EPIC-pn (SW) 

2002/04/08 

2002/04/09 

57.193 

GTO 

effective exposure~71% = 40 ks 

0106260101 

EPIC-MOSl (TI) 

2002/04/08 

2002/04/09 

57.740 

GTO 

no response available yet 

0106260101 

EPIC-MOS2 (SW) 

2002/04/08 

2002/04/09 

57.996 

GTO 


0106260101 

RGSl 

2002/04/08 

2002/04/09 

58.546 

GTO 


0106260101 

RGS2 

2002/04/08 

2002/04/09 

58.546 

GTO 



cal counterpart and its large proper motion of about 
0.33 mas yr“^ as further arguments that RXJ1856 is in¬ 
deed an isolated NS. With data from three HST observa¬ 
tions of RXJ1856, Walter (2001) determined the distance 
to the source to be d = 6llg pc. A re-analysis of the same 
data led to a revised distance to d = 142^3® pc (Kaplan 
et al. 2002). One more HST observation allowed Walter 
& Lattimer (2002) to obtain an improved estimate on the 
distance, d = 117 ± 12 pc. 

None of the previous analysis of the X-ray data yielded 
any significant periodicity. With the ROSAT data Pons et 
al. (2002) derived an upper limit of 6% on the rotational 
modulation of the flux for periods in the range 0.1s - 
20s. BOl, using the data from the first 55-ks Chandra 
observation of RXJ1856, put a limit of 8% in the period 
range 25 ms - 10^ s. Analysing the 505-ks Chandra data 
Ransom et al. (2002) and Drake et al. (2002) obtained 
upper limits of 4.5% (10 ms - 10^ s) and 2.7% (10 ms - 
10^s), respectively. 

In this paper we describe the Chandra and XMM- 
Newton data (§2), the results of spectral (§3) and tim¬ 
ing (§4) analysis, discuss implications on the nature of 
RXJ1856 (§5) and propose alternative models (§6). 


2. Observations and data extraction 

2.1. Chandra LETGS 

RXJ1856 was first observed in March 2000 with the stan¬ 
dard LETGS configuration followed by four further ob¬ 
servations in October 2001 (see Table ||). The LETGS 
spectrum was extracted from the reprocessed level 2.0 
event files. We used the extraction region recommended 
in the Chandra Proposers’ Observatory Guide0 (POG) 
for extracting the source spectrum. For the background, 
large regions above and below the source extraction area 
were selected (see BOl and POG) The measured dis¬ 
persed source count rate is 0.29 ± 0.01 cts s“^ in the 0.15- 
1.00 keV range. The extracted source spectrum and the 

^ http://asc.harvard.edu/udocs/docs/docs.html 


most up-to-date effective area tables (version July 2002)^ 
were used for spectral fits. 

2.2. XMM-Newton 

The XMM-Newton GTO observation took place on April 
8, 2002 (see Table |^). The duration of the exposure for all 
instruments was ~57 ks. The effective exposure for EPIC- 
pn GGD in the Small Window mode is about 70% of this 
giving an effective exposure of ~ 40 ks. The EPIG-pn and 
EPIG-MOS2 countrates in the 0.12-1.2 keV energy range 
are 8.31 ±0.02ctss“^ and 1.35 ±0.01 ctss“^, respectively. 
Two RGS instruments detected the source at the total 
countrate of 0.29 ± 0.03 ctss“^. 


3. Spectral analysis 

BOl found that the RXJ1856 spectrum detected with the 
Chandra LETGS instument in the 55-ks observation is 
best fitted to a blackbody model with kT^ = (63.0 ± 
1.0) eV and x = 4.4 ± 0.5km (d/120pc). A fit to the 
505-ks LETGS spectrum yields very similar, but much 
more constrained values of kT{^y = 63.5 ± 0.2 eV and 
^bbx =4.4±0.1km(d/120pc). The fit is shown in Fig.|| 
(right panels). 

The results of individual fits with blackbody model 
to the spectra from different instruments are given in 
Fig. 1^ and Table^. The Chandra and XMM-Newton data 
yield temperatures well consistent with each other and sig¬ 
nificantly different from that obtained with the ROSAT 
PSPG. The differences in the inferred values of model pa¬ 
rameters are attributed to systematic uncertainties in the 
present calibrations of the X-ray instruments. 


4. Timing analysis 

Details on results from the timing analysis of the Chandra 
data c an be found in BOl, Ransom, Gaensler, & Slane 
(2002) and Drake et al. (2002). The 57-ks observation with 

^ tittp: //asc .harvard.edu/cal/Links/Letg/User 
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Fig. 1. The top panel shows the countrate spectra of RXJ1856 obtained with XMM-Newton and Chandra with the 
best single blackbody model fit to each instrument (the parameters are given in Table |[) . The bottom two panels show 
the ratio between the data and the model for both the CCD detectors and the high resolution grating instruments. 

Table 2. Results of spectral fits to the X-ray data. 


Mission 

Instrument 

Uh 

10^° cm-2 

^-^hh 

eV 

ryoo 

-^bb 

km (d/120 pc) 

/x(0.1 - 1.0 keV) 
10“'^^ erg/cm^/s 

-^bol 

10^^ erg/s 

X/ / d.o.f. 

ROSAT 

PSPC 

1.46±0.20 

56.7±1.0 

7.5±0.5 

1.45 

7.5 

0.9/16 

Chandra 

LETGS 

0.95±0.03 

63.5±0.2 

4.4±0.1 

1.14 

4.1 

1.2/1145 

XMM-Newton 

EPIC-pn 

0.18±0.03 

62.8±0.3 

4.3±0.1 

1.67 

3.7 

2.3/122 

XMM-Newton 

EPIC-MOS2 

0.67±0.02 

62.6±0.4 

4.4±0.1 

1.32 

3.8 

6.1/41 

XMM-Newton 

RGSl-bRGS2 

0.87±0.08 

63.4±0.3 

4.0±0.2 

0.90 

3.3 

1.1/717 


the XMM-Newton EPIC-pn instrument provided much 
larger statistics (370,230 counts extracted from a 30"- 
radius circle centered at the source position) at a 6- 
ms time resolution. The standard Z^-test (Buccheri et 
al. 1983) run in the 10“^ — 50 Hz frequency range re¬ 
vealed a maximum value = 33.3, that translates 

into an upper limit of 1.3% (at a 2a confidence level) on 
variability of the detected radiation assuming a sine-like 
signal. 


5. Discussion 

Pavlov et al. (1996) and Pons et al. (2002) showed that 
the light element (hydrogen and helium) nonmagnetic NS 
atmosphere models can be firmly ruled out because they, 
applied to the X-ray data on RXJ1856, yield (i) too small 
distance estimates {d < 10 pc) and (ii) overpredict the op¬ 
tical flux measured from the source by a factor of ~ 100. 
As first demonstrated by BOl, no acceptable fit can be ob¬ 
tained with iron and standard solar-mixture atmosphere 
models (see Zavlin & Pavlov 2002| for a recent review on 
the NS atmosphere modeling) at any reasonable values 
of gravitational redshift parameter (see BOl for a more 
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Fig. 2. The Chandra LETGS 505 ks X-ray spectrum fitted with a NS solar-mixture atmosphere model (left) and a 
blackbody model (right). From top to bottom the panels show the countrate spectra, the residuals given in sigmas, 
and the photon flux spectra. The model spectra are overplotted in both the top and bottom panels. 


detailed discussion). In case of the solar-mixture compo¬ 
sition, the X-ray data rule out the models with heavy ele¬ 
ment abundances greater than 0.05%. Fig. || (left panels) 
shows an example of a spectral fit with a NS atmosphere 
model where heavy elements provide only 0.2% of the total 
mass density. 

As shown in Zavlin & Pavlov ( 2002| ), sharp features 
in the spectra of NS heavy-element atmospheres, primar¬ 
ily spectral lines, may be smeared out by a fast NS ro¬ 
tation at periods as short as a few milliseconds. However, 
this rotational broadening (the Doppler effect) of the spec¬ 
tral features yet leaves broad-band spectral features in the 
model spectra at energies around most prominent absorp¬ 
tion edges, that make these atmosphere m odels fail to fit 
the X-ray data (Pavlov, Zavlin & Sanwal 2002| ; Braje & 
Romani 2002| ). 

Another possibility would be to fit the data with NS 
atmosphere models for strong magnetic fields. However, 
available magnetized hydrogen models, although not well 
elaborated yet for such rather low temperatures of inter¬ 


est (see Zavlin & Pavlov 2002 for discussion), have the 


same problem as the nonmagnetic case: they overpredict 
the optical flux. On the other hand, the spectra emit¬ 
ted by m agnet ized iron atmospheres (Rajagopal, Romani 
& Miller 1997) show numerous absorption features which 
should be detectable with the modern instruments of high 
energy resolution. We note that smearing of spectral fea¬ 
tures due to line shifts in inhomogeneous magnetic fields 
(for example, varying by a factor of 2 over the NS surface 
for a dipole magnetic held conhguration) may be expected 
to wash away narrow-band features, but anyway, should 
result in broad-band deviations from the blackbody spec¬ 
trum, similar to the case of a fast rotating NS, as discussed 
above. 

We conclude that the “classic” NS atmosphere models 
(with assumption of radiative equilibrium) are unable to 
reproduce the X-ray emission of RXJ1856, which is best 
htted by a simple blackbody model (see Fig. ||, right pan¬ 
els). The possibility of a NS atmosphere which is not in 
radiative equilibrium since its outer layers heated by par¬ 
ticle or photon irradiation has been discussed by Gansicke, 
Braje & Romani (p002). In this case one can produce a 
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Fig. 3. (left) Two-component blackbody model with optical (van Kerkwijk & Kulkarni ^00l| , Pons et al. 2002 ) and 
Chandra LETGS X-ray data, (right) One-component blackbody model that requires a surface emissivity in the X-ray 
band, which is lower than that of a blackbody to fit the X-ray data. In both panels the solid and dotted curves 
represent the absorbed and unabsorbed model spectra, respectively. 


spectrum which is close to that of a blackbody. However, 
this remains a pure speculation until the required source 
of the additional heating is identified. 

An alternative possibility]^ is to assume a condensed 
matter surface — liquid or solid — which might result in 
a virtually featureless Planckian spectrum in the soft X- 
ray band. Such a situation may occur at low temperatures 
{kT < 86 eV) and high magnetic fields {B > 10^^ G) when 
hydrogen, if present on the NS surface, is expected to be 
in t he for m of polyatomic molecules (Lai & Salpeter 1997 ; 
Lai |200l|) . 

Yet another problem arises from the fact that the 
parameters derived from X-rays {kT{^y^ = 63eV and 
^bb.x = 4.4(fi/120pc)km) do not fit the optical data 
obeying the Rayleigh-Jeans law with an intensity about 
a factor of 7 larger than that given by the continuation 
of the blackbody model yielde d by the X-ray data. This 
situation has led Pons et al. (^002 ) to introduce a two- 
component interpretation: the model applied in the X- 
ray band is supplemented with an additional soft black¬ 
body component emitted from about 80% of the NS sur¬ 
face and being responsible for the optical emission (c.f. 
Fig. Ijleft). The requirement that the soft component does 
not contribute in the X-ray band puts an upper limit on 
the blackbody temperature fcTbb.opt < 33.6 eV (a 3cr confi¬ 
dence level), that restricts the stellar radius Rbb,opt > 1^.3 


Originally suggested by G. Pavlov (see BOl for a reference). 


(d/120pc)km (Pavlov, Zavlin & Sanwal 2002). We also 
note that the spectral shape of the optical blackbody data 
puts a lower limit on the temperature of the soft com¬ 
ponent kT^ > 4 eV and an upper limit on its radius 
opt ;£ 46 (d/120 pc) km. 

The non-observed modulation of the X-ray flux im¬ 
poses severe restrictions on the viewing geometry which 
has been discussed quantitatively by Braje & Romani 
( ^0021 ). Using the previous limits on the pulsed fraction 
(<4%) they found ’’the fraction of the sky allowed by 
pulsed fraction constraints” (which translates into a prob¬ 
ability of a given orientation of the rotational and view¬ 
ing axes) of 2-4% for a NS radius of ~ 16 km. The new 
limit established with XMM-Newton on the pulsed frac¬ 
tion reduces the allowed sky fraction to even smaller values 
(around 1%). Alternatively, RXJ1856 could rotate with a 
short period of a few milliseconds (the available X-ray data 
do not provide sufficient time resolution for seaching peri¬ 
odic signals at these time scales). But this seems unlikely 
in view of the NS age of ~ 0.5 Myr implied by its rather 
low surface temperature and the inferred distance from its 
birth place (Walter & Lattimer 2002). 


6. Alternative models 

The absence of periodic variations and spectral features in 
the observed radiation imposes very stringent constraints 
on any model. The simplest way to produce a time con- 
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stant flux would be to assume a uniform temperature 
distribution across the stellar surface. For a temperature 
of kT^ = 63 eV the measured optical spectrum requires 
a blackbody radius i?g^~12.3km. Consequently, the X- 
ray emissivity has to be below that of a blackbody by a 
substantial factor. Using the parameters derived from the 
Chandra data (c.f. Fig. ^ right) we find this factor to be 
about 0.15 (it can be about 0.45 if one adopts the pa¬ 
rameters given by the ROSAT data). This would mean 
that the radiating surface should have a high reflectiv¬ 
ity as may be expected for a condensed matter surface 
(Lenzen & Triimper 1978| ; Brinkmann 1980). In this case 
the spectrum may be represented by a ax x depen¬ 
dence {Bi, is the Planck function), where ax is the ab¬ 
sorption factor (~ [1 — px]j with px being the reflection 
factor), which in the general case will be energy-dependent 
(Brinkmann 1980). We have tested this hypothesis by 
fitting theChandra LETGS spectrum with a Planckian 
Bi, multiplied by an energy dependent absorption factor 
ax = where E is the photon energy. It turns out that 
the best fit yields fd = 1.28 ± 0.30, kT^ = 54 ± 2eV, 
and Tiij = (5.1 ± 0.3) 10^®cm“^. This indicates that at a 
4 cr level we find find deviations from a Planckian spec¬ 
trum which may result from an energy dependent absorp¬ 
tion factor. In this case the radius required from the op¬ 
tical spectrum is ~I3.3a~py^ (d/120pc) km, where 
Oopt < 1 is the absorption factor of the surface in the op¬ 
tical domain. 

In conclusion, the two versions of our one-component 
model yield radii of R°° = 12.3 km and 13.3 km, respec¬ 
tively. These are lower limits as the absorption factor a 
may be smaller than unity in the optical band. We note 
that the observed radii correspond to true NS radii of 
R >9.1km and R >10.3 km (for a NS mass of 1.4 M©), 
respectively, which are consistent with a soft equation of 
state. But, of course, the possible range of parameters al¬ 
lows stiff equations of state as well. 

Two main conditions have to be fulfilled to make this 
model work. Firstly, the NS has to have a condensed mat¬ 
ter surface, which requires a low temperature and a strong 
magnetic field. The former seems to be fulfilled in this 
case: with a temperature of kT = 54 — 63 eV RXJ1856 is 
the coldest one of all detected isolated NSs. But its mag¬ 
netic field is still unknown. The second condition is that 
the condensed matter surface really exhibits the required 
high reflectivity (^ 0.55 — 0.85) in the X-ray domain. It re¬ 
mains to be seen whether a detailed analysis of the optical 
properties of magnetically condensed matter substantiates 
this hypothesis. 


Raumfahrt (BMBF/DLR) and the Max-Planck Society. We 
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comments. 
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Note added in Proof: Turolla, Zane & Drake (2002) have 
submitted a paper to ApJ in which they treat the emissiv¬ 
ity of a condensed matter surface from a theoretical point 
of view, using the method of Brinkman (1980). 
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